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ABSTRACT 

Determining the initial stellar multiplicity is a challenging problem since proto 
stars are faint and deeply embedded at early times; once formed, multiple protostellar 
systems may significantly dynamically evolve before they are optically revealed. In- 
terferometers such as CARMA and ALMA make it possible to probe the scales at 
which turbulent fragmentation occurs in dust continuum emission, potentially con- 
straining early stellar multiplicity. In this Letter we present synthetic observations 
of starless and protostellar cores undergoing fragmentation on scales of a few thou- 
sand AU to produce wide binary systems. We show that interferometric observations 
of starless cores by CARMA should be predominantly featureless at early stages, al- 
though wide protostellar companions should be apparent. The enhanced capabilities 
of ALMA improve the detection of core morphology so that it may be possible to de- 
tect substructure at earlier times. In either case, spatial filtering from inter ferometry 
reduces the observed core substructure and often eradicates traces of existing filamen- 
tary morphology on scales down to 0.025 pc. However, some missing structure may 
be recaptured by combining data from the ALMA full science and Atacama compact 
arrays. 

Key words: stars: formation, stars: starless cores, stars: low-mass 



1 INTRODUCTION 

While the initial mass function of stars is well measured, 
the initial stellar mul tiplicity remains largely unconstrained 
(Duchene et al. 2007|). Determining multiplicity among field 
stars well after the formation process has ended is itself a 
challenging problem: close binaries are often unresolvable, 
gravitational boundedness is difficult to confi rm, and back- 
grou nd stars are ever-present interlopers ([Raghavan et al.l 
2010). During the formation and subsequent evolution, dy- 
namical interactions influence initial companion separations 
and decrease multiplicity. 

Observing the initial protostellar multiplicity in situ is 
challenging for similar reasons with the added complica- 
tion that young protostars are dim and heavily obscured by 
dust and gas. However, recent interferometric instruments 
such as the Combined Array for Research in Millimeter- 
wave Astronomy (CARMA) and the Atacama Large Mil- 
limeter/submillimeter Array (ALMA) are revolutionizing 
the ability to probe the earliest stages of cores on scales 
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of a few arcseconds or less. This presents an opportunity 
to observe core substructure, protostellar disks, and young 
companions. 

There are two main theories of binary star formation. In 
the disk fragmentation scenario, massive protostellar accre- 
tion disks become Toom re unstable and fragment into one or 
more close companions (| Adams et al . 1989; B onnell &; Bate! 
Il994l ). In the turbulent fragmentation scenario, turbulent 
perturbations within a single prestellar core or filament in- 
dividuall y collapse to form separate stars with wide sepa- 
ratio ns ([Goodwin et al.l I2QQ4I : [Fisherl 1200^ : iGoodwin et al.l 
l2007h . Numerical and analytic arguments indicate that disk 
insta bility should be quite common in high-mass star forma- 
tion (jKratter fe Matznerl [20061 : iKratter et al.l l2QQ8h . How- 
ever radiation feedback significantly reduces the fragmenta- 
tion of disks around low-mas s stars, resultin g in few low- 
mass multiple st ar systems (|Cai et al.l 120081 : lOffner et al.l 
120091 : lBateN2009h . The multiple systems that do form are 
the result of turbulent frag mentation an d have initial sepa- 
rations of -1000 AU (jOffner et al.ll2010h . 

The enhanced capabilities and sub-arcsecond resolution 
of ALMA will make it possible to test both scenarios of bi- 
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Figure 1. Single dish observation of 1.1 mm flux in mJy per pixel (0.13") of a binary protostellar system formed via turbulent 
fragmentation, where the system is placed at 250 pc. Left: Perfect resolution simulation data. Center: Image convolved with a 5" 
beam. Right: Image convolved with a 31" beam. The image extent is 80" across at a distance of 250 pc. Crosses mark the protostar 
positions. Contours indicate 30%, 50%, 70% and 90% of the image maximum. 



nary formation and map protostellar accretion disks down 
to AU size scales. Several authors have used interferometric 
obser vations to look for young companions of Class objects 
(e.g.. lLoonev et~aT1l20od : IjQrgensen et alJl2007l : iMaurv et all 
2010). However, they arri ve at somewhat different conclu- 
sions: lLoonev et al.l (|2000h find that al l their target e d em - 
bedded objects have companions, while IMaurv et al.l ((20101) 
find only one tent ative companion i n a sa mple of five Class 
sources. Likewise, Jprg ensen et al.l (|2007h identified only one 
candidate previously unknown companion. Using continuum 
emission, it is possible in principle to detect companions 
forming as a result of tu rbulent fragm e ntatio n at even ear- 
lier times. For example, ISchnee et al.l (|2010l ) use CARMA 
to observe 3mm continuum emission from 11 starless cores 
in the Perseus molecular cloud at 5" resolution. They found 
that the cores had no conspicuous substructure. This could 
be due to several possibilities. The cores may be too young 
or may never go on to form an individual star much less a 
wide multiple system. Turbulent fragmentation may be very 
uncommon. Or turbulent fragmentation may be ongoing in 
the cores, but beyond the observable limits of CARMA. We 
investigate the third possibility here. 

If cores, such as those investigated by ISchnee et al.l 

(|2010h . actually contain young protostars, it might be pos- 
sible to observe fragmentation at shorter wavelengths. In- 
deed, a number of very cold cores thought t o be starless 
have since been found to contain protostars ([Enoch et al.l 
l20ld : iDunham et ai1l201ll : IPineda et aDl201lh . However, in 
the case of L1451, the presence of a protostar, potentially 
still in the "first core" gas stage, was ultimately identified 
by outflow activity rather than by thermal emission. This 
reinforces the point that source identification is challeng- 
ing, and it requires instruments with high sensitivity and 
resolution to characterize core structure and identify young 
companions. 

In this Letter, we use the CAS A software package 
to synthetically observe binaries forming due to turbulent 
fragmentation in the numerical simulations of Offner et al. 
(2009). By following the evolution of such pairs beginning 



in the prestellar core stage, we can make predictions about 
the feasibility of observing such fragmentation in dust con- 
tinuum and constraining stellar multiplicity at the earliest 
stage of star formation. 



2 NUMERICAL SIMULATIONS AND 
METHODS 

The molecular clou d simulations we obs erve in this letter are 
those presented in I Offner et al.l (|2009h . Since Offner et al. 
(2009) fully describe the calculations, we include only a 
brief overview here. The calculations are performed using 
the ORION adaptive mesh refinement code including driven 
large-scale turbulence, self-gravity, and flux-limited diffusion 
radiative transfer. Forming stars are modeled by sink par- 
ticles with a sub-grid stellar evolution model. These are in- 
serted in regions of the flow that exceed the maximum grid 
resolution. In practice, since protostellar winds are not in- 
cluded, these particles give an upper limit on the stellar 
mass. The cloud domain size is 0.65 pc across with a mini- 
mum cell size of 32 AU. 

In our analysis, we focus on several typical systems out 
of the approximately 10 close pairs that form over the course 
of the simulation. Since the systems form self-consistently 
from the turbulent gas, the fragmentation history and stellar 
masses are not predetermined. 

We use the "simdata" task in the Common Astronomy 
Software Applications (CASA) packag^B to produce syn- 
thetic interferometric observations of the starless and proto- 
stellar systems. Using simdata, the model cloud was placed 
at the RA and Dec of the Perseus molecular cloud. For the 
CASA input, we require maps of the cores in units of flux. 
Converting between simulated column density and contin- 
uum flux is straight-forward for optically thin emission. As- 
suming a constant dust temperature, Td, the flux at a given 
frequency, can related to the column density by: 

1 http://casa.nrao.edu 
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Figure 2. Comparison of simulations and CARMA synthetic observation for fragmenting cores at three different times. Left: Simulation 
data with perfect resolution in 1.1mm flux (Log mjy beam -1 ). Color scale is the same as in Figure 1. Right: CARMA observation of 
each core at 3mm emission including Gaussian as noise assuming a distance of 250 pc in mJy beam -1 . Contours indicate 30%, 50%, 
70% and 90% of the maximum flux. The CARMA beamsize is indicated by the hatched ovals. 



S v = EB„(T D )/c A n 6 , (1) 

where B u (Td) is the Planck function evaluated at the dust 
temperature, k\ is the dust opacity, Qb is the beam solid 
angle and £ is t he gas column density per pixel (e.g., 
lEnoch et al.l [2006). The opacities for 1.1 mm and 3 mm 
are aci.i 0.0114 cm 2 g" 1 and k 3 = 0.00169 cm 2 g" 1 
(|Ossenkopf fe Hennind 11994). The simulations assume that 
the dust and gas are well-coupled, whi ch is a reasonable ap- 
proxi mation for densities > 10 4 cm -3 ([Goldsmith &; Langerl 
1 19781 ). Thus, we adopt Td — 10 K, the simulation gas tem- 
perature, w hich is also a lower limit for the dust temperature 
in Perseus (| Schnee et al J [20091 ). Once protostars form, they 
heat their environment to mean temperatures of 15-20 K 
and result in higher emission. 

Producing synthetic single dish observations of the sim- 
ulations is a straight- forward application of equation [1] Fi- 
nite resolution can then be imposed by convolving the flux 
map with a circular Gaussian beam. 



3 SYNTHETIC OBSERVATIONS 

For the purpose of comparison, we first present synthetic 
single dish observations of the fragmenting cores. We then 
produce interferometric observations mimicing the specifi- 
cations of CARMA and ALMA. Finally, we investigate the 
influence of noise and distance on structure detection. 

3.1 Single Dish Observations 

Observations that map out entire molecular clouds and iden- 
tify star-forming cores often use single dish continuum data. 



However, with beam resolutions of tens of arcseconds, re- 
solving core substructure in even nearby clouds is impossi- 
ble. Figure [1] shows a protostellar system forming in a dense 
core with 0" (i.e., perfect), 5" and 31" resolution. The two 
peaks are distinct with 5" resolution, although the fainter 
peak may not be apparent depending upon sensitivity and 
noise levels. At 31", a resolution comparable to SCUBA (850 
fim) and Bolocam (1.1 mm), only a single peak is apparent. 

Interferometers, such as CARMA, are able to achieve 
this 5" scale, but can only probe a fixed window of scales. 
Since larger scale information is resolved out, target fluxes 
may also be reduced by 90% or more relative to single dish 
observations. Figure [1] illustrates that cores forming wide 
protobinary companions may be observable with current in- 
terferometer technology. 



3.2 CARMA Observations 

We produce synthetic o bservations to comp are with the 
CARMA observations of Schnee et al.1 (|20ld ). In order to 
make a truly similar comparison, we reproduce their observ- 
ing procedures and conditions as closely as possible. Our ob- 
servations have a 2.8 GHz bandwidth centered at 102 GHz, 
and we adopt the CARMA D-array configuration and an 
integration time of 100 seconds with a total time of 8 hours. 
Unless otherwise stated, we assume the systems are 250 pc 
away: the distance of the Perseus molecular cloud. The pixel 
size of the simulated map is 0.13" , and the synthesized beam 
is 5" . We add synthetic Gaussian noise with as = 0.7 m Jy 
beam -1 , comparable to that of the Schne e et al. | J2010l ) ob- 
servations. 

Figure [2] shows synthetic CARMA observations of three 
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Figure 3. Comparison of early science (ES) and full science (FS) 
ALMA configurations for Core 1 (left) at t\ and Core 3 (right) 
at £3 assuming a distance of 250 pc. Color scale indicates flux 
in mJy beam -1 . The top two panels are observed with the early 
science (ES) extended configuration; the bottom two panels are 
observed in one example full science (FS) configuration. White 
contours indicate 30%, 50%, 70% and 90% of the maximum flux. 
Orange contours indicate 10% of the maximum flux. The ALMA 
beamsize is indicated by the hatched ovals. 

separate cores undergoing turbulent fragmentation at three 
different times. The simulation data seen in 1.1 mm emission 
is plotted for comparison. At the first time all three cores 
are starless. However, each shows some evidence of collapse 
and fragmentation at perfect resolution. The second row il- 
lustrates that fragmentation proceeds on the order of 10-20 
kyr, which is a relatively short time compared to the typical 
~100 kyr dynamical time of cores. If the core exists in a 
quasi - steady state for some time before undergoing collapse 
(e.g. jBroderick et al.ll2008f h then the likelihood of catching 
any particular starless core in the act of fragmenting may be 
small. Such short timescales for o bserving close companion s 
are consistent with those found bv lStamatellos et al.l (|201l[ ), 
who synthetically observed massive fragmenting disks. The 
probability is reduced further since not all cores (may) ex- 
perience fragmentation on ~ 1000 AU scales. 

At the earliest times, interferometric observations do 
not clearly show fragmentation. The second fragment in 
Core 1 becomes more apparent over time, but is nearly in- 
visible in the starless phase. However, the low level detec- 
tions do appear similar in size and separation to some of 
the lower fl ux contours appare nt in maps of Perbo45 and 
Perbo58 bv lSchnee et al.l (|2010|). In contrast, all of the fila- 
mentary structure in the starless Cores 2 and 3 is resolved 
ou t. These resul t s are consistent with the core observations 
of ISchnee et al ] (|2010h . who find that most starless cores 
that have bright 1.1mm emission in single dish maps are 
undetected in 3mm interferometric maps. This also confirms 
that additional filamentary structure may be removed by the 
interferometric technique. 

The synthetic observations show that at later times the 
protostars and companions become brighter. Due to the fac- 
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Figure 4. Synthetic observation combining data from ALMA FS 
with ALMA compact array for Core 1 (top) at t\ and Core 3 
(bottom) at £3. Contours indicate fluxes with values 10% (or- 
ange), 30%, 50%, 70% and 90% of the maximum. The effective 
beamsize is indicated by the hatched oval. The observations as- 
sume 100s integration time and two hours of observing with each 
configuration. 



tor of 10 difference in opacity, the over-densities will be sig- 
nificantly brighter in 1.1mm, which will increase the signal 
to noise. This suggests that widely separated protostellar 
companio ns should be relativel y apparent at high resolu- 
tion fe.g. jMerrill fe EnocrJl201oh . However, the length of the 
window in which fragmentation occurs is still problematic if 
secondaries migrate to shorter separations or are unbound 
on short timescales. 



3.3 ALMA Predictions 

ALMA, which will have 66 reconfigurable antennas, sub- 
arcsecond angular resolution and sensitivity to wavelengths 
from 3mm to 300 /xm, will significantly expand observational 
capabilities. Although currently only partially completed, 
the ALMA early science program allows reduced observa- 
tions with 16 12m operational antennas. Figure [3] shows 
simulated observations of Core 1 and Core 3 comparing the 
ALMA early science Cycle extended configuration and the 
ALMA full science configuration 10. We adopt configuration 
10 in the ALMA full science library in simdata since it is a 
representative intermediate configuration. 

Each observation assumes a 100 second integration time 
and a total time of 2 hours. We add thermal noise to the 
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Figure 5. Left: Simulated 3mm CARMA observation of Core 1 at t\ observed at a distance of 250 pc with 2<rs noise. Center: Simulated 
3mm CARMA observation of Core 1 observed at a distance of 500 pc with la noise. Right: Simulated 3mm full science ALMA observation 
of Core 1 observed at a distance of 500 pc. Color scale shows flux in mjy beam -1 . Contours indicate fluxes with values 30%, 50%, 70% 
and 90% of the maximum. The beamsize is indicated by the hatched ovals. 



simulated ALMA images using simdata, assuming that the 
precipitable water vapor (PWV) during the observations is 
2.8mm (appropriate for Band 3 ALMA observations) and 
a bandwidth of 8 GHz. ALMA maps were cleaned until 
a threshold of 3 times the rms thermal noise in the map 
was reached. The pixel size, as for the CARMA simulations, 
is 0.13", and the early science and full science synthesized 
beam sizes are ^3" and 1.5", respectively. 

Figure [3] illustrates that ALMA better resolves the pri- 
mary peak morphology. A secondary peak in Core 1 is 
marginally shown by contours with 10% of the maximum 
flux. However, the filamentary structure of Core 3 is not 
visible. The early science configuration sensitivity to sub- 
structure is otherwise similar to that of CARMA with four 
times the observing time. The ALMA full science configu- 
ration detects the Core 3 filament with some suggestion of 
further fragmentation. This highlights the difficulty of imag- 
ing structure and fragmentation even with ALMA's superior 
resolution. 

It is possible to increase the range of recoverable spatial 
scales by combining data from antenna configurations with 
different baselines. The size of the ALMA 12-meter anten- 
nas prohibit them from being placed within 15 meters of one 
another, which limits the maximum resolvable spatial scale. 
To permit a greater range of observations, the completed 
ALMA site will include a second smaller array, the Ata- 
cama Compact Array (ACA), comprised of four 12-meter 
and 12 7-meter antennas. Here, we combine the higher spa- 
tial resolution synthetic data of the intermediate FS main 
array configuration with synthetic data from the 7-meter 
ACA antennas. Figure [4] shows a simulated observation of 
two cores where the visibilities from FS and ACA have been 
added and then deconvolved using the CASA CLEAN sub- 
routine. The appearance of filamentary structure and core 
substructure is significantly improved compared to the sin- 
gle configuration data in Figures [2] and [3] 



3.4 Noise and Resolution Limitations 

Object distance and observation sensitivity are both critical 
to mapping core structure. Figure [5] shows Core 1 with in- 
creased noise and source distance. At our fiducial noise and 
resolution (see Figures [2] and [3]) both increased noise and 
distance eliminate any detection with CARMA. At a dis- 
tance of 500 pc, the secondary peak is also undetected with 
the ALMA full science configuration. However, the primary 
peak is reliably detected. This implies that fragmentation at 
later times, such as that illustrated in Core 1 at 50 kyr and 
Core 3 at ^ 9 kyr, would be readily observable. 



4 CONCLUSIONS 

In this Letter we produce synthetic observations of frag- 
menting starless and protostellar cores. We show that inter- 
ferometric observations of starless cores by CARMA should 
be predominantly featureless at early stages. In fact, struc- 
ture may be apparent only within a short period, ~ lOkyr, 
of the formation of a protostar. This may account for some 
of the apparent l ack o f substructure in starless cores noted 
bv lSchnee et al.1 (|2010h . We find that wide protostellar com- 
panions with separations of ~ 1000AU should be detectable. 
The confidence of the secondary detection depends upon the 
source age, resolution, and signal-to- noise, which may par - 
tially explain the diffe ring r esults of lLoonev et al.l (|2000h , 
IjQrgensen et al.1 (|2007T ). and iMaurv et al.1 (feOloj TALMA's 
enhanced capabilities improve the detection of core mor- 
phology, so that it may be possible to detect substructure 
at earlier times. Filamentary structure is more difficult to 
detect than peakiness, and interferometry, especially at high 
resolution, significantly reduces the presence and extent of 
filaments. However, we find that it is possible to recover 
missing structure by combining ALMA Full Science data 
with data from the Atacama Compact Array. 
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